Diagnosis, prescription and prognosis of a Bell-state filter by quantum process 

tomography. 
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Using a Hong-Ou-Mandel interferometer, we apply the techniques of quantum process tomography 
to characterize errors and decoherence in a prototypical two-photon operation, a singlet-state filter. 
The quantum process tomography results indicate a large asymmetry in the process and also the 
required operation to correct for this asymmetry. Finally, we quantify errors and decoherence of the 
filtering operation after this modification. 

PACS numbers: 42.50-p, 03.67.Mn, 03.67.Pp 



Quantum computation promises exponential speedup 
in the solution of difficult problems such as factoring large 
numbers and simulating quantum systems 0, 0j- In a 
quantum computer single- and multiple-qubit operations 
drive the system through a sequence of highly entangled 
states before the result is finally measured. A quantum 
computation is vulnerable to errors and to environmen- 
tal decoherence, which destroys the entanglement. Char- 
acterization of quantum operations including errors and 
decoherence is a pressing issue for quantum information 
processing and is possible by the technique of quan- 
tum process tomography (QPT) |4L|5|]. QPT has been 
demonstrated for single qubits la, |j| and for mixed en- 
sembles of two-qubit systems in NMR Here we 
present QPT of an entanglement-generating two-qubit 
operation, the partitioning of photons by a beamsplitter 
in a Hong-Ou-Mandel (HOM) interferometer. Our char- 
acterization reveals large imperfections in the process and 
indicates the appropriate remedy. Finally, we extend the 
QPT results to predict the accuracy of the process, once 
repairs are carried out. 

Multi-qubit operations on photons, once thought to 
require very large optical nonlinearities, can now be per- 
formed with linear optical elements such as wave-plates 
and beamsplitters coupled with the highly nonlinear pro- 
cess of photodetection. This idea is exploited in schemes 
for linear optics quantum computation HoL 111! Il2l and 
to generate multi-photon entangled states [TJ, |hj| . The 
schemes are probabilistic and employ post- selection: the 
photodetection signals indicate when the correct opera- 
tion has taken place. The HOM effect plays a central 
role in these proposals, and itself is a prototypical exam- 
ple of a post-selected multi-qubit operation; it generates 
correlations and entanglement without optical nonlinear- 
ities. The HOM effect has been used toproduce entan- 
gled states for Bell inequality tests 0, and to make 
probabilistic Bell state measurements for quantum tele 
portation and entanglement swapping 
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they are in some way distinguishable [19j. We use pho- 
ton pairs which are indistinguishable in wavelength, spa- 
tial mode and arrival time at the beamsplitter, leaving 
only the polarization to (possibly) distinguish them. By 
detecting photons leaving from different output ports, 
we post-select an entangled polarization state. Ideally, 
the process acts as a filter for the Bell singlet state 
^~ = (HV — VH)/\/2, in which the photons have or- 
thogonal polarizations in any basis. In any real appa- 
ratus this process will include errors and decoherence. 
Using the techniques of QPT, we determine how the po- 
larization state, more specifically the 4x4 density matrix 
p which describes an arbitrary two-photon mixed state, 
changes in passing the beamsplitter. In general, p will 
change as p( m > — > p( out > = £(p' m )), where £ is the super- 
operator, a linear mapping from input density matrices to 
output density matrices. The superoperator completely 
characterizes the effect on the system, including coherent 
evolutions, decohering interactions with the environment, 
and loss. 
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In the HOM effect, two photons meeting at a 50/50 
beamsplitter can leave by different output ports only if 



FIG. 1: Schematic of experimental setup. BBO, /3-barium 
borate crystals, H half-wave plate, Q quarter wave plate, 
BS beamsplitter, PBS polarizing beamsplitter, SPCM single- 
photon counting module, TRR translatable retro-reflector. 



We use a HOM interferometer constructed to produce 
arbitrary input polarizations and detect arbitrary output 
polarizations. The experimental setup is shown schcmat- 
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ically in Fig. 1. A 7 mW beam of 351.1 nm light from 
an argon-ion laser illuminates a pair of 0.6 mm thick 
/3-barium borate crystals, cut for degenerate downcon- 
version at a half-opening angle of 3.3° . Pairs of down- 
conversion photons at 702.2 nm emerge from the crystals 
vertically polarized. This initial polarization state can be 
rotated into any input product state by the state prepara- 
tion half- and quarter- wave plates immediately before the 
central beamsplitter. The downconversion beams meet 
the beamsplitter at 45° incidence. The beamsplitter it- 
self [13] consists of a multi-layer dielectric coating on a 
glass substrate, with an anti-reflection coated back face. 

Polarization analyzers consisting of a quarter- and a 
half-waveplate before a polarizing beamsplitter are used 
to select an arbitrary product state. Photons which 
pass the analyzers are detected by single-photon counting 
modules and individual and coincidence detection rates 
are registered on a computer. Downconversion beams 
were aligned to overlap both spatially and temporally on 
the beamsplitter, giving a HOM dip visibility of 90 ± 5% 
for both horizontal and vertical input polarizations. The 
process tomography measurements described below were 
performed at the center of this dip. 

We prepare 16 linearly independent input states 
{Pi™* 1 } an d measure the corresponding outputs {pf wt ^}- 
The inputs [2l| are the pure states = where 
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{tpi, . . . ,ipie} = 



{HH, HV,VV,VH, 
RH, RVDV, DH, 
DR, DD, RD, HD, 
VD,VL,HL,RL} 



(1) 



and the polarizations are horizontal H, vertical V , diag- 
onal D = (H + V)/V%, right circular R = (H - iV)/\?2 
and left circular L = (H +iV) / y/2. A single output pf wt ^ 
can be found by making projective measurements onto 
the sixteen states {ipj} ■ The coincidence rates for these 

measurements are Ra = 



i?oTr[p[ out) |^-)(^|], where R 
is the constant rate of downconversion at the crystals. 
Note that we use non-normalized output density matri- 
ces, i.e. Tr[/?(°"*'] < 1, because photon pairs can be lost in 
the process. Absorption and scattering losses are small, 
but post-selection necessarily removes a significant frac- 
tion of the pairs for most input states. 

The measured coincidence rates Rij are shown in Fig. 
2. As expected for a filter, the output has similar po- 
larization characteristics for all inputs, but not all are 
equally transmitted, e.g., HH and VV are blocked. 
A typical output density matrix, reconstructed using 
maximum-likelihood estimation [2l| is shown in Fig. 3. 
The large coherence between HV and VH indicates that 
this is an entangled state, with a concurrence 
of C = 0.89. The HOM effect is acting as an entangled- 
state filter, but the selected state is clearly not which 
has a real density matrix and negative off-diagonal ele- 




HH VV RH DV DR RD VD HL 
HV VH RV DH DD HD VL RL 



FIG. 2: Coincidence rates. Brightness indicates the count rate 
observed in a given two-photon polarisation state (horizontal 
axis) for a given input polarisation state (vertical axis). 
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FIG. 3: Output density matrix (normalized) for an input 
state of HV. Left graph shows Re[p'°"*'], right graph shows 
Im[p (out) ]. 



ments. 

We can understand this behaviour through the su- 
peroperator £. For clarity, we work in the Bell state 
basis where * ± = (HV ± VH)/^/2 

and $ ± = (HH ± VV)/y/2. We use a matrix rep- 
resentation for £: The density matrix is written as a 
real 16-dimensional vector p made from the independent 
coefficients of the non-normalized density matrix, i.e., 
P = (pn, ■ ■ ■ ,P44,Re[pi 2 ],Im[pi 2 ],Re[pi3], . . .,lm[p 34 ]) T . 
The superoperator is represented by a matrix M which 
acts as 



(2) 



In principle, M could be found from this equation by a 
simple inversion, since we measured Rij for a basis set 
{p[ m ^}- This procedure is sensitive to small errors and 
can produce a non-physical M, i.e., one which predicts 
non-physical (mathematically, non-positive semidefinite) 
p(out) ^ Instead, we reconstruct M by maximum- 
likelihood estimation within the space of completely posi- 
tive superoperators, i.e. operators that map physical den- 
sity matrices to physical density matrices (see e.g. Su- 
darshan [24| for the mathematical conditions on the map- 
ping between density matrices). The reconstructed M is 
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FIG. 4: Reconstructed superoperators for the post-selected 
HOM process, a) superoperator as measured, b) predicted 
superoperator after repair. The matrix M is shown, input 
density matrix elements at bottom, output elements at left, 
where the density matrices are represented in vector form (see 
text). Horizontal stripes on the vertical bars indicate the the 
best estimated value and the statistical uncertainties. 

shown in Fig. 4 a), with error estimates from an en- 
semble of simulated datasets Poisson distributed around 
the measured data. This matrix is "normalized" to give 
Tr[p out ] = 1/4 when the input is a completely mixed 
state. 

We verify the accuracy of the reconstructed superop- 
erator using the the input states LL and RR, which were 
not used in the reconstruction process. These states are 
used as input, both to equation J2J an d in the HOM 
interferometer. In both cases, prediction and the exper- 
iment result (again by maximum likelihood reconstruc- 
tion) agree with fidelity of 97%. 

The superoperator M bears little resemblance to an 
ideal singlet-state filter, for which My = S^iSj^. Clearly 
the process is not performing the intended filtering op- 
eration. In fact, it is nearly a projection onto a different 



maximally entangled state |25|. Written as a canonical 
Kraus operator sum j2(| [2^] , the superoperator allows us 
to find this state directly. In the sum £{p) = J2i KipKj , 
the leading operator K\ is very nearly a projector onto 
the state = (HV - exp[i4>]V H) / with 4> = 0.84 tt. 
This immediately suggests a way to "fix" the non-ideal 
beamsplitter. Adding a birefringent phase shifter which 
takes VH — > exp[— i<fi]VH before the beamsplitter and 
the reverse operation afterward would give (nearly) a 
single-state filter. 

Finally, we can predict the behaviour of this "fixed" 
operation. The corresponding matrix M is shown in Fig. 
4b). The large (1,1) element indicates the filtering opera- 
tion and the smaller nonzero elements contribute to deco- 
herence and other errors. These errors presumably arise 
from imperfect overlap of the downconverted beams and 
residual imperfections in the beamsplitter. They do not 
appear to have a simple form, but we can gain some in- 
sight from some simple measures, calculated using the su- 
peroperator. An unpolarized input (a completely mixed 
state) gives rise to an output that is 84% 'F - , or an av- 
erage polarization ratio (intensity of versus average 
intensity of the other three Bell states) of 16:1. This same 
output is entangled, with a concurrence of C = 0.70, suf- 
ficient for a Bell inequality violation. Quantifying purity 
with the linear entropy Sl, which ranges from for a 
pure state to 1 for a completely mixed state, the process 
purifies the mixed state from Sl = 1 to Sl = 0.37. We 
can also ask how well the repaired filter maintains the co- 
herence of an input. The pure input , F _ is passed 75% of 
the time and emerges largely pure, with Sl = 0.13. The 
state \E' + is 13% passed with Sl — 0.51 and <F are on 
average 6% passed with low purity Sl = 0.88. Of course, 
different applications for a singlet-state filter will have 
different requirements and different figures of merit. The 
superoperator we have found using QPT is more general, 
a complete characterization suitable for evaluating any 
proposed use. It is also, we have seen, a useful diagnostic 
and predictive tool. 
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